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Abstract. 

I review the status of observational determinations of central masses in 
nearby galactic nuclei. Results from a variety of techniques are summarized, 
including ground-based and space-based optical spectroscopy, radio VLBI 
measurements of luminous water vapor masers, and variability monitoring 
studies of active galactic nuclei. I will also discuss recent X-ray observations 
that indicate relativistic motions arising from the accretion disks of active 
nuclei. The existing evidence suggests that supermassive black holes are an 
integral component of galactic structure, at least in elliptical and bulge- 
dominated galaxies. The black hole mass appears to be correlated with the 
mass of the spheroidal component of the host galaxy. This finding may have 
important implications for many astrophysical issues. 



1. Motivation 

The discovery of quasars in the early 1960's quickly spurred the idea that 
these amazingly powerful sources derive their energy from accretion of mat- 
ter onto a compact, extremely massive object, most likely a supermassive 
black hole (SMBH; Zel'dovich & Novikov 1964; Salpeter 1964; Lynden-Bell 
1969) with M 10^ - 10^ M©. Since then this model has provided a 
highly useful framework for the study of quasars, or more generally, of the 
active galactic nucleus (AGN) phenomenon (Rees 1984; Blandford & Rees 
1992). Yet, despite its success, there is little empirical basis for believing 
that this model is correct. As pointed out by Kormendy & Richstone (1995, 
hereafter KR), our confidence that SMBHs must power AGNs largely rests 
on the implausibility of alternative explanations. To be sure, a number of 
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characteristics of AGNs indicate that the central engine must be tiny and 
that relativistic motions are present. These include rapid X-ray variability, 
VLBI radio cores, and supcrluminal motion. However, solid evidence for 
the existence of SMBHs in the centers of galaxies has, until quite recently, 
been lacking. 

As demonstrated by Soltan (1982), simple considerations of the quasar 
number counts and standard assumptions about the efficiency of energy 
generation by accretion allows one to estimate the mean mass density 
of SMBHs in the universe. The updated analysis of Chokshi &; Turner 
(1992) finds p, 2 x IO^Cq J M0 Mpc^"^ for a radiative efficiency of 
e = O.leo.i- Comparison of p, with the B-hand galaxy luminosity den- 
sity of 1.4xl0^/i Lq Mpc~^ (Lin et al. 1996), where the Hubble constant 
Hq = WOh km s~^ Mpc~^, implies an average SMBH mass per unit stel- 
lar luminosity of ~1.4x 10~^eQ Mq/Lq. A typical bright galaxy with 
L% « 10^°/i"^ Lq potentially harbors a SMBH with a mass > WeQ\h~^ 
Mq. These very general arguments lead one to conclude that "dead" quasars 
ought to be lurking in the centers of many nearby luminous galaxies. 

The hunt for SMBHs has been frustrated by two principal limitations. 
The more obvious of these can be easily appreciated by nothing that the 
"sphere of influence" of the hole extends to ~ GM,/a'^ (Peebles 1972; 
Bahcall Wolf 1976), where G is the gravitational constant and a is 
the velocity dispersion of the stars in the bulge, or, for a distance of D, 
~l"(M./2 X 10^ Mq)(o-/200 km s~^)~'^{D/5 Mpc). Typical ground-based 
observations are therefore severely hampered by atmospheric seeing, and 
only the heftiest dark masses in the closest galaxies can be detected. The 
situation in the last few years has improved dramatically with the advent 
of the Hubble Space Telescope (HST) and radio VLBI techniques. The more 
subtle complication involves the actual modeling of the stellar kinematics 
data, and in this area much progress has also been made recently as well. 

Here I will highlight some of the observational efforts during the past 
two decades in searching for SMBHs, concentrating on the recent advances. 
Since this contribution is the only one that discusses nuclear BHs aside 
from that in the Milky Way (Ozcrnoy, these proceedings) and in NGC 
4258 (Miyoshi, these proceedings), I will attempt to be as comprehensive 
as possible, although no claim to completeness is made, as this is a vast 
subject and progress is being made at a dizzying pace. To fill in the gaps, 
I refer the reader to several other recent review papers, each of which has 
a slightly different emphasis (KR; Rees 1998; Richstone 1998; Ford et al. 
1998; van der Marel 1998). 
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2. Early Clues from Photometry 

The prospect of finding massive BHs in globular clusters motivated much 
early effort to investigate the distribution of stars resulting from the adi- 
abatic growth of a BH in a preexisting stellar system. The central density 
deviates strongly from that of an isothermal core and instead follows a 
cuspy profile p{r) oc r~^/^ (Young 1980) or steeper if two-body relaxation 
(Peebles 1972; Bahcall & Wolf 1976) or different initial density profiles 
(Quinlan, Hernquist, & Sigurdsson 1995) are taken into account. The dis- 
covery that the centers of some giant elliptical galaxies obey this prediction 
generated much enthusiasm for the existence of SMBHs. In the well-known 
case of M87 (Young et al. 1978), Lauer et al. (1992) have since shown that 
the central cusp persists to the limit of the resolution of the HST (C'l). 

However, as emphasized by Kormendy (1993; see also KR), photomet- 
ric signatures alone do not uniquely predict the presence of a SMBH. The 
cores of most galaxies are now known to be nonisothermal. And more- 
over, contrary to naive expectations, galaxy cores with high central surface 
brightnesses and small core radii, far from being the ones most likely to 
host SMBHs, are in fact least expected to do so. This apparently con- 
tradictory statement can be most easily understood by considering the so- 
called fundamental-plane relations for the spheroidal component of galaxies 
(Faber et al. 1987; Bender, Burstein & Faber 1992). More luminous, more 
massive galaxies tend to have more massive central BHs (§ 7), but they 
also have larger, more diffuse cores. Indeed, high-resolution photometric 
studies of early-type galaxies (Nieto et al. 1991; Crane et al. 1993; Jaffe et 
al. 1994; Lauer et al. 1995) find that the central surface brightness profiles 
cither continue to rise toward the center as I{r) cx r~'^, with 7 ^ 0.5- 
1.0 (the "power-law" galaxies) or they fiatten at some characteristic radius 
to a shallower slope of 7 « 0.0-0.3 (the "core" galaxies). The power-law 
galaxies are invariably lower luminosity, lower mass systems compared to 
those with distinct cores. 

In summary, photometric signatures alone cannot be used as reliable 
indicators for the presence of SMBHs. Instead, we must turn to the more 
arduous task of obtaining kinematic measurements. 

3. Methods Based on Stellar Kinematics 

Contrary to the ambiguity of light profiles, the Keplerian rise in the ve- 
locity dispersion toward the center, a{r) oc r~^/^, is a robust prediction 
for a wide variety of dynamical models containing a central massive dark 
object (MDO; Quinlan et al. 1995). Sargent et al. (1978) noticed that the 
innermost velocities of M87 were consistent with such a prediction, and, 
assuming an isotropic velocity distribution, they inferred that the center 
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of this galaxy contained a dark mass of ~5xl0 M.q, presumably in the 
form of a SMBH. The central rise in a{r), unfortunately, can be insidiously 
mimiccd by an anisotropic velocity distribution, and therefore an MDO is 
not required by the data for this object (Duncan & Wheeler 1980; Binney & 
Mamon 1982; Richstone &; Tremaine 1985; Dressier & Richstone 1990; van 
der Marel 1994a). This degeneracy presents a serious difficulty for many 
mass determinations based on stellar kinematic data. An extensive and lu- 
cid discussion of this vast subject was presented by KR, and many of the 
details will not be repeated here. Nonetheless, an abbreviated synopsis is 
needed to motivate the topic. 

Following the notation of KR, the radial variation in mass can be ex- 
pressed by the first velocity moment of the collisionless Boltzman equation, 

M(r) ^Y^^^\_^_ ^) 
^ ' G ^ G [ dlnr dlnr ^ a^' ^ . 

where V is the rotational velocity, is the radial and ag and the 
azimuthal components of the velocity dispersion, and v is the density of the 
tracer population. In practice, several simplifying assumptions are adopted: 
(1) the mass distribution is spherically symmetric; (2) the mean rotation is 
circular; and (3) v is proportional to the luminosity density, or, equivalently, 
that M/L does not vary with radius. 

A brief scrutiny of the above equation indicates that the effects of ve- 
locity anisotropy can have a large and complicated effect on the derivation 
of M(r) because the terms inside the bracket significantly affect the a'jr/G 
term. If Gr > ere and cr^ > cr^, each of the last two terms will be neg- 
ative and can be as large as -1. The central brightness distributions of 
the spheroidal component of most galaxies typically have — (dln^-Zdlnr) 
-1-1.1 for luminous, nonrotating systems and t +2 for low to intermediate- 
luminosity systems (e.g., Faber et al. 1997). Since — (din cr^/dln r) < +1, 
it is apparent that, under suitable conditions, all four terms can largely 
cancel one another. As emphasized by KR, all else being equal, smaller, 
lower luminosity galaxies such as M32 potentially yield more secure mass 
determinations than massive galaxies like M87 because less luminous sys- 
tems tend to have (1) steeper central light profiles, (2) a greater degree of 
rotational support, and (3) less anisotropy. 

The principles behind the stellar kinematics analysis are conceptually 
straightforward but in practice technically challenging. Given the set of ob- 
served quantities I{r), V{r), and cr{r), the goal is to derive a range intrinsic 
values for these quantities after accounting for projection and the blurring 
effects of seeing. Much of the machinery for these tasks has been devel- 
oped and extensively discussed by Kormendy (1988a, b) and Dressier & 
Richstone (1988). The sensitivity of the results to the effects of anisotropy 
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are examined through maximum-entropy dynamical models (Richstone &; 
Tremaine 1984, 1988) to see whether conclusions regarding the presence 
of MDOs can be obviated by a suitable exploration of parameter space. 
Perhaps the most serious limitation of these maximum-entropy models is 
that they do not properly take flattening into account. 

The last several years have seen a resurged interest in improving the 
techniques of analyzing stellar kinematics data. In the context of SMBH 
searches, Gerhard (1993), van der Marel et al. (1994a, b), Dehnen (1995), 
among others, have stressed the importance of utilizing the full informa- 
tion contained in the velocity profile or line-of-sight velocity distribution 
(LOSVD) of the absorption lines, which are normally treated only as Gaus- 
sians. A system with significant rotation, for instance, can leave a measur- 
able skewness on the LOSVD, while various degrees of anisotropy would 
imprint symmetric deviations from a Gaussian line shape. Neglecting these 
subtleties can lead to systematic errors in the measurement of V{r), but 
in the cases best studied so far these effects do not seem to have been se- 
vere (KR). Furthermore, the line profile should develop weak, high- velocity 
wings if a SMBH is present (van der Marel 1994b), although the currently 
available data do not yet have the requisite quality to exploit this tool. 

Yet another advance has focused on the development of dynamical mod- 
els with two-integral phase-space distribution functions, f{E,Lz), E being 
the total energy and the angular momentum in the symmetry axis (van 
der Marel et al. 1994b; Qian et al. 1995; Dehnen 1995). Such models are 
properly flattened, and they generate predictions for the LOSVDs; on the 
other hand, it is not clear whether imposing a special dynamical structure 
is too restrictive. This limitation will be eliminated by fully general, ax- 
isymmetric three-integral models (van der Marel et al. 1998; Cretton et al. 
1998; Gebhardt et al. 1998). 

There are currently 10 galaxies with published MDO measurements 
determined from stellar kinematical data (Table 1). Of these, only three 
(M81: Bower et al. 1996; NGC 3379: Gebhardt et al. 1998; NGC 4342: van 
den Bosch 1998) come solely from HST data; the remaining ones, although 
many by now confirmed with HST^ were initially discovered from high- 
quality ground-based observations (sec KR for a detailed account of each 
object). Kormendy and collaborators, in particular, making use of the excel- 
lent seeing conditions and instrumentation on the CFHT, continue to make 
progress in this area. Two new MDOs have been reported recently based on 
CFHT data: the low-luminosity elliptical galaxy NGC 4486B has Mmdo = 
6x10^ M© (Kormendy et al. 1997b), and NGC 3377, another close cousin, 
has Mmdo = 2.3x10^ M© (Kormendy et al. 1998). This demonstrates the 
important fact that even in the HST era ground-based observations con- 
tinue to play an important role in SMBH searches. 
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Figure 1. Stellar kinematic data for NGC 3115 compared with various dynamical models 
(see Kormendy et al. 1996). The left panel shows the best ground-based data, and the 
right panel the same data with new HST data superposed for comparison. Both V and 
a rise much more steeply with radius in the new HST data. 



The new observations with HST^ thus far all acquired using the Faint 
Objeet Spectrograph (FOS), provide an important contribution by increas- 
ing the angular resolution by about a factor of 5 compared to the best 
ground-based data available. In all cases studied (NGC 3115: Kormendy 
et al. 1996; NGC 4594: Kormendy et al. 1997a; M32: van der Marel et al. 
1997; M31: Ford et al. 1998), the velocity dispersions continue to rise to- 
ward smaller r and the maximum rotational velocity has generally increased 
(Fig. 1). In the case of NGC 3115, the FOS spectra are of sufficient quality 
to reveal wings in the LOSVD that extend up to ~1200 km s"^ (Kor- 
mendy et al. 1996). The HST data thus considerably bolster the case for 
a MDO in these objects. The improvement in angular resolution addition- 
ally strengthens our confidence that the MDOs might indeed be SMBHs. 
A reduction of the size scale by a factor of 5 increases the central density 
by more than two orders of magnitude. Although in general this is still not 
enough to rule out alternative explanations for the dark mass (§ 6), it is 
clearly a step in the right direction. 

I conclude this section with a few remarks on the dark mass in the 
Galactic Center (see Ozernoy in these proceedings for more details), which, 
in my view, is now the most compelling case of a SMBH in any galactic 
nucleus. Prom analysis of an extensive set of near-IR radial velocities of in- 
dividual stars, coupled with additional measurements from the literature. 
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Figure 2. Enclosed mass versus radius for the Galactic Center derived from stellar radial 
velocities and proper motions (from Genzel et al. 1997). The points labeled with "G" 
come from gas kinematics. The thick dashed curve is a model for the stellar cluster with 
a total mass density of 4x10® M© pc~^; the solid curve denotes the sum of this cluster 
and a point mass of 2.61x10® M©; and thin dashed curve is the sum of the stellar cluster 
and an additional dark cluster with a central density of 2.2x10^^ M© pc~^. 

Genzel et al. (1996; see also Krabbe et al. 1995) found a highly statisti- 
cally significant rise in the radial velocity dispersion between 5 and 0.1 
pc from the dynamical center. Assuming an isotropic velocity distribution, 
the observations require a dark mass of ~3xl0^ Mq within r = 0.1 pc 
and M/Lk > 100; the dark mass must have a density in excess of 10^ 
Mq pc~^, which argues strongly for it being a SMBH. These conclusions, 
and a suspicion nearly three decades old (Lynden-Bell & Rees 1971), have 
finally been vindicated by recent measurements of stellar proper motions 
within the central 1 pc region using high-resolution K-hand astrometric 
maps (Eckart & Genzel 1996, 1997; Genzel et al. 1997; Ghez et al. 1998). 
The main results are the following: (1) the stellar radial velocities agree 
with the proper motions, which implies that on average the velocities are 
close to isotropic; (2) the combined velocities imply a dark mass (Fig. 2) 
within 0.006 pc of 2.61 x 10^ M© (Genzel et al. 1997 quote a statistical error 
of ±0.15 and a combined statistical and systematic error of ±0.35); (3) the 
density, therefore, has an astonishingly high value of >2xl0"'^^ Mq pc~"^, 
which leaves almost no room to escape the conclusion that the dark mass 
must be in the form of a SMBH (§ 6). The presence of a large mass is 
also supported by the detection of several stars, within 0.01 pc from the 
central radio source Sgr A*, moving at speeds in excess of 1000 km s~^. 
Prom the velocities of the fast-moving stars and the near stationarity of 
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Sgr A*, Genzel et al. further use equipartition arguments to constrain the 
mass of the radio core itself (> 10^ M©), which, when combined with the 
exceedingly small upper limit for its size (r < 4xl0~^ pc), would imply a 
density of >3xW^^ M© pc'^. 

4. Methods Based on Gas Kinematics 

Unlike the situation for stars, gas kinematics are much easier to interpret 
if the gas participates in Keplerian rotation in a disklike configuration. But 

there are two caveats to remember. First, gas can be easily perturbed by 
nongravitational forces (shocks, radiation pressure, winds, magnetic fields, 
etc.). Indeed, in the case of the Galactic Center, it was precisely this reason 
that its central mass, which had been estimated for some time using gas 
velocities (Lacy ct al. 1980), could not be accepted with fiill confidence 
prior to the measurement of the stellar kinematics. Second, there is no a 
priori reason that the gas should be in djmamical equilibrium, and therefore 
one must verify empirically that the velocity field indeed is Keplerian. The 
optically-emitting ionized gas in the central regions of some spirals show 
significant noncircular motions (e.g., Fillmore, Boroson, 8z Dressier 1986). 
NGC 4594 is a striking example. Kormendy et al. (1997a) showed that 
the emission-line rotation curve near the center falls substantially below 
the circular velocities of the stars, and hence the gas kinematics cannot be 
used to determine the central mass. 

4.1. OPTICAL EMISSION LINES 

The sharpened resolution of the refurbished HST has revealed many ex- 
amples of nuclear disks of dust and ionized gas (Fig. 3) . The nuclear disks 
typically have diameters ~100-300 pc, with the minor axis often aligned 
along the direction of the radio jet, if present. Some examples include the 
elUptical galaxies NGC 4261 (Jaffe et al. 1993), M87 (Ford et al. 1994), 
NGC 5322 (Carollo et al. 1997), and NGC 315 (Ho et al. 1998), and the 
early-type spiral M81 (Devereux, Ford, Sz Jacoby 1997). I will highlight here 
only three cases; Table 1 gives a complete list of objects and references. 

The first object for which the nuclear gas disk was used to determine 
the central mass was M87. Harms et al. (1994) used the FOS to obtain 
spectra of several positions of the disk and measured a velocity difference 
of ±1000 km s~^ at a radius of 0'.'25 (18 pc) on either side of the nucleus. 
Adopting an inclination angle of 42° determined photometrically by Ford 
et al. (1994), the velocities were consistent with Keplerian motions about 
a central mass of (2.4±0.7) x 10^ Mq. Since the implied M/Ly ~ 500, 
Harms et al. concluded that the central mass is dark, most likely in the 
form of a SMBH. The case for a SMBH in M87 has been considerably 
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strengthened through the recent reobservation with HST by Macchetto et 
al. (1997), who used the long-sht mode of the Faint Object Camera to 
obtain higher quahty spectra extending to r = O'.'OS (3.5 pc). The velocities 
in the inner few tenths of an arcsecond are well fitted by a model of a thin 
disk in Keplerian rotation (Fig. 4), although the inclination angle is not 
well constrained (47°-65°). The rotation curve at larger radii falls below 
the Keplerian curve, possibly indicating a warp in the disk (Macchetto et 
al. 1997) or substantial perturbations due to spiral shocks (Chakrabarti 
1995). For i = 52°, Mmdo = (3.2±0.9) x 10^ M©, and M/Ly > 110. If, 
instead, a Plummer potential is assumed, the distributed dark mass can 
have a core radius no larger than ~5 pc. So, in either density ~10 

M.Q pc^^ is implied. 

The mildly active nucleus of NGC 4261 contains a rotating disk of dust 
and ionized gas as well (Ferrarese, Ford, & Jaffe 1996); like M87, the disk 
is slightly warped and shows traces of weak spiral structure. Although the 
FOS data for this object are rather noisy, they indicate that the gas largely 
undergoes circular motions. The mass interior to r = 15 pc is Mmdo = 
(4.9±1.0)xl0® M0, and M/Ly has an exceptionally high value of 2x10^. 
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Figure 4- Optical emission- line rotation curve for the nuclear disk in M87. The two 
curves in the upper panel correspond to Keplerian thin disk models, and the bottom 
panel shows the residuals for one of the models (see Macchetto et al. 1997). 



The installation of the imaging spectrograph STIS in 1997 at long last 
gives HST an efficient means to obtain spatially resolved spectra of the 
central regions of galaxies. Much progress in the field is anticipated in the 
near future. A taste of what might be expected can be seen in the early- 
release observations of M84 by Bower et al. (1998; Fig. 5). M84 is almost 
a twin of M87 in terms of luminosity, and its central dark mass (1.5x10^ 
M0), too, is similar. 

Lastly, I mention an interesting, unconventional case. The radio galaxy 
Arp 102B belongs to a minority of AGNs that display so-called double- 
peaked broad emission lines. Several ideas have been proposed for the pe- 
culiar line profiles in this class of objects, but the favored explanation is 
that the lines originate from a rclativistic accretion disk (Eracleous et al. 
1997). During the course of a long-term optical monitoring of Arp 102B, the 
intensity ratio of the two peaks of the Ha line displayed sinusoidal varia- 
tions with a period of 2.2 years for several years (Newman et al. 1997). The 
periodic signal was interpreted as arising from a "hot spot" in the accretion 
disk. By modeling the line profile from the epochs when the hot spot was 
quiescent, one can estimate the radius and inclination angle of the spot's 
orbit, and, combined with its period, the enclosed mass. The mass within 
r = 0.005 pc turns out to be 2.2x10^ Mq, consistent with a moderately 
luminous {Mb ~ -20 mag) elliptical (see § 7). 
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Figure 5. Optical cmission-linc rotation curve for the nuclear disk in M84 obtained 
with STIS on HST (Bower at al. 1998). The kinematics can be well fitted with a thin 
Keplerian disk model, which is plotted as open circles connected by the solid line. 

4.2. RADIO SPECTROSCOPY OF WATER MASERS 

Luminous 22-GHz emission from extragalactic water masers are preferen- 
tially detected in galaxies with active nuclei, where physical conditions, 
possibly realized in a circumnuclear disk (Claussen & Lo 1986), evidently 
favor this form of maser emission. With the detection in NGC 4258 of 
high-velocity features offset from the systemic velocity by ~ ±900 km s"^ 
(Nakai, Inoue, &; Miyoshi 1993), Watson & Wallin (1994) already surmised 
that the maser spectrum of this Seyfert galaxy can be interpreted as aris- 
ing from a thin Keplerian disk rapidly rotating around a mass of ~10^ M©. 
But the solid proof of this picture came from the high-resolution (A^ = 
0'.'0006x0'.'0003; Aw = 0.2 km s'^) VLBA observations of Miyoshi et al. 
(1995) who demonstrated that the maser spots trace a thin (<0.003 pc), 
nearly edge-on annulus with an inner radius of 0.13 pc and an outer radius 
of 0.26 pc. The systemic features lie on the near side of the disk along the 
line-of-sight to the center (Fig. 6); the high- velocity features delineate the 
edges of the disk on either side and follow a Keplerian rotation curve to 
very high accuracy (< 1%). The implied binding mass within 0.13 pc is 
3.6x10^ Mq, which corresponds to a density of >4xl0^ Mq pc~^. In fact, 
one can place a tighter constraint on the density. The maximum deviation 
of the velocities from a Keplerian rotation curve limits the extent of the 
central mass to r ^ 0.012 pc (Maoz 1995), from which follows that the 
density must be >5xl0^^ Mq pc~^. 

Two, possibly three, additional AGNs have H2O megamasers suitable 
for tracing the central potential. The spectrum of the maser source in the 
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Figure 6. Water maser emission in NGC 4258 (Miyoshi et al. 1995). Top: spatial distri- 
bution of the maser features; bottom: rotation curve. Adapted from Greenhill (1997). 



Seyfert nucleus of NGC 1068 also exhibits satellite features (±300 km s^^) 
offset from the systemic velocity (Greenhill et al. 1996). The redshifted and 
blueshifted emission again lie on a roughly linear, 2-parsec feature passing 
through the systemic emission (Greenhill 1998). The rotation curve in this 
instance is sub-Keplerian, possibly because the disk has nonnegligible mass, 
and the derived mass (1.7x10"^ Mq within a radius of 0.65 pc) is less certain. 

The maser in the nucleus of NGC 4945 shows a position and velocity 
distribution reminiscent of NGC 4258 as well: high-velocity features sym- 
metrically straddle the systemic emission. Greenhill, Moran, & Herrnstein 

(1997) interpret the data, which in this case is considerably less accurate 
because of its location in the southern sky, in terms of an edge-on disk 
model and derive a central mass of 1.4x10^ M0 within r = 0.3 pc. This 
result is quite surprising because, as an Scd spiral, NGC 4945 is expected 
to be essentially bulgeless. If the dark mass in its center is truly in the form 
of a SMBH, then SMBHs evidently can form without a bulge. 

The H2O megamaser source in NGC 3079 is potentially useful for mass 
determination. Here, however, the complex spatial distribution of the emis- 
sion regions and the large intrinsic widths of the lines complicate the anal- 
ysis, and the interpretation of the data may not be unique. Trotter et al. 

(1998) tentatively assign a central mass of 1x10^ Mq to this galaxy. 
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4.3. DETERMINING CENTRAL MASSES OF ACTIVE GALACTIC NUCLEI 

I mention one other method for determining masses in the central regions 
of galaxies, specifically in AGNs. Although AGNs largely provide the mo- 
tivation for searching for SMBHs, ironically it is precisely in these objects 
that conventional techniques used to measure masses fail. The bright con- 
tinuum emission of the active nucleus nearly always completely overpowers 
the stellar absorption lines near the center, and in many cases the narrow 
emission lines are significantly affected by nongravitational forces. 

An approach taken in the past attempts to utilize the broad [(1-few) x 10^ 
km s^^] emission lines that are thought to arise from the so-called broad-line 
region (BLR) , a tiny, dense region much less than a parsec from the central 
source. Assuming that the line widths trace gravity, the mass follows from 
iIV^'TB'Lk/G^ where r] ~ 1-3 depending on the kinematic model adopted. 
The BLR radius has traditionally been estimated from photoionization ar- 
guments (e.g., Dibai 1981; Wandel & Yahil 1985; Wandel & Mushotzky 
1986; Padovani, Burg, & Edelson 1990), but recent variability studies indi- 
cate that the BLR is much more compact than previously thought (Netzer 
& Peterson 1997). 

The continuum output from AGNs typically varies on timescales rang- 
ing from days to months in the UV and optical bands. Because the emis- 
sion lines arc predominantly photoionized by the central continuum, they 
vary in response to the changes in the continuum, but with a time delay 
(lag) that corresponds to the light-travel distance between the continuum 
source and the line-emitting gas. "Reverberation mapping" (Blandford & 
McKee 1982), therefore, in principle allows one to estimate the luminosity- 
weighted radius of the BLR, although in practice the complex geometry 
and ionization structure of the BLR complicate the interpretation of the 
"sizes" derived by this method (see Netzer & Peterson 1997 for a recent 
review) . 

If the widths of the broad emission lines reflect bound gravitational 
motions, as seems to be the case in most well-studied objects (Netzer & 
Peterson 1997; but sec Krolik 1997), then, adopting a reasonable kinematic 
model (e.g., randomly moving clouds), the virial mass can be estimated 
from v^t-slr/G. If, instead, the clouds are infalling, as has been claimed in 
some cases, the mass will be smaller by a factor of 2. One of the uncertainties 
in the application of this simple formalism lies in the choice of v. What 
is appropriate? One reasonable choice might he v = (\/3/2)FWHM, the 
full width at half-maximum of a representative broad line. Yet another 
ambiguity is which line to use, since not all broad emission lines have the 
same widths. Ultraviolet or high-ionization lines, for instance, generally 
have broader profiles than optical or low-ionization lines. For the purposes 



TABLE 1 

MASSES DETERMINED FROM STELLAR AND GAS KINEMATICS 
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Type 


(Mpc) 


(mag) 
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(1) 
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(4) 


(5) 




(6) 


(7) 
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Milky Way 


Sbc 
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006 
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0.30 


6.6x10* 


>100 


g 


Q 


M32 


E2 


0.72 


8.72 
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0.23 
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G 
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14 
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NGC 4945 
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0.05 
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4.0 


8.32 






-14.51'" 


<4xl0° 


11.0 






S 


25 


Arp 102B 
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REFERENCES. — (1) Kormendy & Richstone 1995; (2) Kormendy &; McClure 1993; (3) Filippenko & Ho 1998; (4) Simien & de Vaucouleurs 1986; (5) Genzel et al. 1997; (6) Ford et al. 
1998; (7) van der Marel et al. 1997; (8) Bower et al. 1996; (9) Bower et al. 1998; (10) Macchetto et al. 1997; (11) Jones et al. 1996; (12) Greenhill 1998; (13) Trotter et al. 1998; (14) 
Kormendy et al. 1996; (15) Kormendy et al. 1998; (16) Gebhart et al. 1998; (17) Miyoshi et al. 1995; (18) Ferrarese, Ford, & Jaffe 1996; (19) van den Bosh 1998; (20) Kormendy et al. 
1997b; (21) Kormendy et al. 1997a; (22) Greenhill, Moran, & Herrnstein 1997; (23) Ferrarese. Ford. ^ Jaffe 1008: (2-1) van der Marel 1008; (25) Maiolino et al. 1998; (26) Newman et al. 
1997. 

Note. — Cols.: (l) Galaxy name; (2) Hubble type from de Vaucouleurs et al. 1991; (3) distance from Tully 1988 or derived by assuming Hq — 75 km s~ Mpc~ ; (4) total apparent 
blue magnitude corrected for Galactic and internal extinction, from de Vaucouleurs et al. 1991; (5) ratio of bulge to total luminosity based on the photometric studies given in the adjacent 
column; (6) absolute blue magnitude of the bulge component; (7) mass of the MDO, taken from the study given in the last column; (8) radius contaning Mmdo; (9) density derived from 
radius and Affjoo; (lO) mass-to-light ratio (in the V band) of the MDO; (ll) method used: "G" — gas kinematics; "S" — stellar kinematics; "V" — variability. 

^Upper limit on size based on arguments by Maoz 1995. 
Bulge luminosity estimated assuming M/Lb — 6, where M — bagRe/G is the mass of the bulge, uq is the central velocity dispersion, and is the effective radius. From Maiolino 
et al. (1998), ctq 80 km s~^ and R^ 80 pc. 
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of this exercise, I simply chose the hne for which the most data exist (H/9) 
in order to obtain as large a sample as possible. 

TABLE 2 

MASSES DETERMINED FROM REVERBERATION MAPPING" 



Galaxy 


Hubble 


D 


B°. 


B/T 


MB(bul) 


rsLR 


FWHM(H,3) 


Mmdo 




Type 


(Mpc; 


1 (mag) 




(mag) 


(ltd) 


( km s-') 
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(6) 


(7) 


(8) 
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13.64 


0.86 


-21.61 


39 
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Fairall 9 


S? 
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13.50 


0.08 


-20.13 


23 


4200 
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SBb 


89 


13.32 


0.20 
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18 
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1.0x10* 
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Pair? 
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20 
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0.12 
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13.66 


0.47 


-20.63 


21 
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S? 
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14.30 


0.50 
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16 


4100 


3.8x10' 


NGC 3227 


SABa 


20.6 


11.18 


0.52 


-19.68 


17 


3900 


3.8x10' 


NGC 3516 


SBO: 
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12.14 


0.61 


-20.27 


7 


4760 
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NGC 3783 


SBa 


38.5 
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0.33 


-19.68 
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NGC 4151 


SABab 


20.3 
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0.36 


-19.72 
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NGC 4593 


SBb 


39.5 


11.43 


0.48 


-20.76 


4 
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NGC 5548 


SO/a 


67.0 


12.81 


0.47 


-20.50 


19 


5610 
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NGC 7469 


SABa 


65.2 


12.64 


0.40 


-20.44 


5 


3388 


9.1x10" 


Note. — Cols. 


: (1) Galaxy 


name; 


(2) Hubble type 


from de Vaucouleurs 


et al. 1991; (3) distance from TuUy 1988 or 



derived by assuming ^fo — 75 km Mpc~^; (4) total apparent blue magnitude corrected for Galactic and internal 
extinction, from de Vaucouleurs et al. 1991; (5) ratio of bulge to total luminosity (host galaxy plus AGN); Mrk 110 is 
too disturbed to yield a reliable photometric decomposition; (6) absolute blue magnitude of the bulge component; (7) lag 
(in light days) between the continuum and the H,/3 hght curves; (8) FWHM of the broad emission hne; (9) virial mass 
derived from tblr and FWHM(H/3). 

^All Scyfcrt 1 nuclei studied to date with reverberation mapping of the hne. Optical monitoring data exist for 3C 
390.3, but it was not included because the profiles of its Balmer lines are very complicated. References for the data entries 
have been omitted for brevity. 



Table 2 lists the derived masses for the 17 Seyfert 1 galaxies that have 
been monitored extensively in the optical; eight of the objects appear in the 
compilation of Peterson et al. (1998). Since the masses of MDOs derived 
from gas and stellar kinematics show a loose correlation with the bulge or 
spheroidal luminosity of the host galaxies (§ 7; Fig. 8a), we can ask whether 
those derived from reverberation mapping follow such a correlation. I have 
estimated the i?-band luminosities of the bulges of the Seyferts based on 
published surface photometry of the host galaxies (taking care to exclude 
the contribution of the AGN itself, which often can be significant). Figure 
86 indicates that, at a fixed bulge luminosity, the masses from reverberation 
mapping are systematically lower than the masses obtained using conven- 
tional techniques, on average by about a factor of 5. It is encouraging that 
this admittedly crude method of mass estimation is not too far off the mark. 
Notably, the scatter of Mmdo at a fixed luminosity is quite comparable in 
the two samples, and the constant offset suggests that one of the under- 
lying assumptions in the mass estimate is incorrect. Since the line width 
affects the mass quadratically, it is conceivable that some measure of the 
line profile other than the FWHM is more appropriate. 
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Rest energy (keV) 

Figure 7. The Fe Kq line in the composite spectrum of Seyfert 1 nuclei (Nandra et al. 
1997). The solid line is a fit to the line profile using two Gaussians, a narrow component 
centered at 6.4 keV and a much broader, redshifted component. 



5. Indirect, but Tantalizing Evidence 

Lastly, one additional piece of evidence, although it does not give a direct 
measure of the central mass, cannot go unmentioncd — namely the recent 
detection in AGNs of the broad iron Ka line at 6.4 keV. This line has been 
known for some time to be a common feature in the hard X-ray spectra of 
AGNs, and it is thought to arise from fluorescence of the X-ray continuum 
off of cold material, presumably associated with the accretion disk around 
the SMBH (e.g.. Pounds et al. 1990). The spectral resolution of the existing 
data, however, was insufficient to test the predicted line profile (Fabian et 
al. 1989). The ASCA satellite provided the much-awaited tell-tale signs in 
the deep exposure of the Seyfert 1 galaxy MCG-6-30-15 (Tanaka et al. 
1995): the Fe Ka line exhibits Doppler motions that approach relativistic 
speeds (~100,000 km s~^ or 0.3c!) as well as an asymmetric red wing 
consistent with gravitational redshift. The best-fitting disk has an inner 
radius of 6 Schwarzschild radii. The relativistic Fe Ka line, now seen in 
a large number of sources (Nandra et al. 1997; Fig. 7), provides arguably 
the most compelling evidence to date for the existence of SMBHs. Other 
mechanisms for generating the line profile are possible, but implausible 
(Fabian et al. 1995). Detailed modeling of the line asymmetry has even the 
potential to measure the spin of the hole, but this is still very much a goal of 
the future given the current data quality and uncertainties in the modeling 
itself (e.g., Reynolds & Begelman 1997; Rybicki & Bromley 1998). 

6. Are the Massive Dark Objects Really Black Holes? 

Thus far we have rigorously shown only that many galaxies contain cen- 
tral MDOs, not that the dark masses must be in the form of SMBHs. 
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Direct proof of the existence of SMBHs would require the detection of rel- 
ativistic motions emanating from the vicinity of the Schwarzschild radius, 
Rs = 2GM,/c^ 10-^(M./10S M©) pc. Even for our neighbor M31, 
Rs subtends 3xl0~^ arcseconds, and the Galactic Center only a factor of 
2 larger. We are clearly still far from being able to achieve the requisite 
angular resolution and in the meantime must rely on indirect arguments. 

One approach seeks to identify some observational feature that might be 

taken as a fingerprint of the event horizon or of physical processes uniquely 
associated with the environment of a BH. One such "signature" might be 
the broad Fe Kck line discussed in § 5; another is the high-energy power- 
law tail observed in some AGNs and Galactic BH candidates (Titarchuk & 
Zannias 1998). And yet a third possibility is the advection of matter into 
the event horizon (Menou, Quataert, & Narayan 1998). 

A different strategy appeals to the dynamical stability of the probable 
alternative sources of the dark mass (Goodman k. Lee 1989; Richstone, 
Bower, & Dressier 1990; van der Marel et al. 1997; Maoz 1998). The ab- 
sence of strong radial gradients in the stellar population, as measured by 
variations in color or spectral indices, implies that the large increase in 
M/L toward the center cannot be attributed to a cluster of ordinary stars. 
On the other hand, the undcrluminous mass could, in principle, be a clus- 
ter of stellar remnants (white dwarfs, neutron stars, and stellar-size BHs) 
or perhaps even substellar objects (planets and brown dwarfs). To rule out 
these possibilities, however exotic they might seem, one must show that the 
clusters cannot have survived over the age of the galaxy, and hence finding 
them would be highly improbable. 

As most recently discussed by Maoz (1998), the two main processes 
that determine the lifetime of a star cluster are evaporation, whereby stars 
escape the cluster as a result of multiple weak gravitational scatterings, 
and physical collisions among the stars themselves. Exactly which domi- 
nates depends on the composition and size of the cluster, and its maximum 
possible lifetime can be computed for any given mass and density. Maoz 
(1998) shows that in two galaxies, namely the Milky Way and NGC 4258, 
the density of the dark mass is so high (> 10^^ M© pc~^) that it cannot 
possibly be in the form of a stable cluster of stellar or substellar remnants: 
their maximum ages [~(l-few)xlO^ yr] are much less than the ages of the 
galaxies. The only remaining constituents allowed appear to be subsolar- 
mass BHs and elementary particles. This constitutes very strong evidence 
that the MDOs — at least in two cases — are most likely SMBHs. In the 
following discussion, I will adopt the simplifying viewpoint that all MDOs 
are SMBHs, bearing in mind that at the current resolution limit we cannot 
yet disprove the dark-cluster hypothesis for the majority of the objects. 
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7. The Black-Hole Mass/Bulge Mass Relation 

Does M, depend at all on other properties of the host galaxies? A much- 
discussed possibility is that M, scales with the mass of the spheroidal com- 
ponent of the host (Kormendy 1993; KR; Faber et al. 1997; Magorrian et al. 
1998; Richstone 1998; Ford et al. 1998; van der Marel 1998). The significance 
of the scatter in the correlation, or whether any correlation exists at all, is 
not yet certain. It is somewhat disconcerting that different authors plotting 
the same objects do not always arrive at the same conclusion. The discrep- 
ancies can often be traced to different assumptions about distances, source 
of bulge-to-disk decomposition, and even apparent magnitudes adopted for 
the host galaxies (e.g., extinction is not always corrected). The set of host 
galaxy parameters I adopt is compiled in Table 1. 

Figure 8 a illustrates that there indeed appears to be a trend of M, 
increasing with bulge mass (luminosity). It is encouraging to note that 
the central masses derived from gas and stellar kinematics do not show any 
obvious systematic offsets relative to one another. No obvious differentiation 
by Hubble type is evident either. As has been noted by others, the scatter 
of M, at a given luminosity is considerable, at least a factor of 10, perhaps 
up to 100. The scatter may have been exacerbated slightly by four possibly 
anomalous points. NGC 4486B is a companion to M87, and it appears to 
have been tidally truncated; its original luminosity was probably higher. 
On the other hand, the bulge luminosity of NGC 4945 could very well have 
been overestimated. Its bulge-to-disk ratio was found using the relation of 
Simicn & dc Vaucouleurs (1986), which may be inappropriate for a galaxy 
of such late Hubble type (Scd). Finally, the masses of M81 and NGC 3079 
are quite uncertain and probably have been underestimated. 

The trend is much more significant when five upper limits are included. 
NGC 205, a dwarf elliptical companion of M31, contains a blue, compact 
nucleus with characteristics resembling an intermediate-age globular clus- 
ter. Its core radius, determined from HST photometry, combined with a 
ground-based measurement of its velocity dispersion yields an upper limit 
of 9x10^ Mq for any dark mass (Jones et al. 1996). The bulgcless, late-type 
(Scd) spiral M33 also has a stringent upper limit on its central mass. Its 
nuclear cluster is extremely tiny (core radius ^ 0.39 pc), and its central 
velocity dispersion is 21 km s^^; Kormendy k. McClure (1993) put an up- 
per limit of M. < 5x10^ M0. NGC 4395 in some ways resembles M33, but 
it is even more extreme. The nucleus is optically classified as a type 1.8 
Seyfert (broad Ha and H/? present), emits a largely nonstellar featureless 
continuum that extends into the UV (Filippenko, Ho, & Sargent 1993), 
and displays variable soft X-ray emission and a compact flat-spectrum ra- 
dio core (Moran et al. 1998). These properties alone would be unremarkable 
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Figure 8. (a) Log Af. versus log Ls(bulgc) for the objects listed in Table 1. The typical 
uncertainty of M, is probably about a factor of 2. Open symbols denote points that 
may have an exceptionally large uncertainty in either of the two variables (see text). 
Masses derived from stellar kinematics are plotted as circles, those from gas kinematics 
as squares, the unconventional case of Arp 102B as a triangle, and five upper limits as 
crosses. Lines of constant mass to luminosity ratio axe also shown. (6) Same as in (a), 
but for the Seyfert galaxies listed in Table 2 (except Mrk 110). 
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were it not for the fact that the nucleus has an absolute blue magnitude of 
only -9.8 and lives in a Magellanic spiral 2.6 Mpc away! Filippenko &; Ho 
(1998) detected the Ca II infrared triplet lines in absorption from echelle 
spectra taken with the Keck telescope, from which they were able to es- 
timate the strength of the stellar component contributing to the nuclear 
light {Mb = -7.3 mag) and the central stellar velocity dispersion {a ~ 30 
km s~^). Combining the velocity dispersion with a cluster size (r ^ 0.7 pc) 
obtained from HST images, Filippenko &: Ho limit the central mass to ^ 
8x10^ Mq. The Circinus galaxy is thought to house a Seyfert nucleus, and 
if it contains a SMBH, its mass within r ~ 10 pc has been constrained to 
be < 4x10^ Mq (Maiohno et al. 1998). The last upper limit shown in the 
figure pertains to the globular cluster M15; following KR, I adopt an upper 
limit of M, = 1x10^ M0. 

However, before reading too much into this diagram, we should ask 
whether the apparent correlation might arise from selection effects. The 
absence of points on the upper left-hand corner is probably real; there 
is nothing preventing us from detecting a massive BH in a small galaxy. 
Yet, we should be cautious, because very few low-mass galaxies have been 
studied so far, most of the effort having been focused on luminous, early- 
type systems. On the other hand, the empty region on the lower right- 
hand corner could be an artifact. Small masses are difficult to detect at 
large distances, and most luminous galaxies are far away. So the apparent 
correlation could be an upper envelope. Future observations are needed to 
settle this issue. 

The median value of M,/LB(bul) for the 20 detected objects is 0.012, 
which translates into a mass ratio of 0.002 for M/Lb ~ 6 typical for 
old stellar populations (van der Marel 1991). That is, on average about 
0.2% of the bulge mass is locked up in the form of a SMBH. Magorrian 
et al. (1998) constructed axisymmctric f{E,Lz) models for a sample of 32 
early-type (mostly E and SO) galaxies having both HST photometry and 
ground-based stellar kinematics data, and they concluded that the data are 
consistent with nearly all of the galaxies having SMBHs. The 29 detected 
objects have a median M,/Mbui ~ 0.005, higher than found here. However, 
as Magorrian et al. realize, the assumption of a two-integral distribution 
function may have caused them to overestimate M, (cf. van der Marel 
1998). Interestingly, quasars possibly also obey a similar M,-Mbui relation. 
McLeod (1998) finds that, for the most luminous quasars, there exists a 
minimum host luminosity that increases with nuclear power. Assuming that 
the quasar luminosities correspond to energy generation at the Eddington 
rate, M,/Mbui is again ~0.002 (McLeod 1998). 

With regard to the dead quasar prediction discussed in § 1, recall that 
we expect to find on average a 10^ Mq BH for every Lb ^ 10^^ Lq galaxy. 
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Figure 9. (o) Detection rate and (6) number distribution of AGNs as a function of 
Hubble type in the spectroscopic survey of Ho et al. (1995, 1997). "Type 1" AGNs (those 
with broad Ha) are shown separately from the total population (types 1 and 2). 



or M,/Ls(bul) « 3.3xl0~^ Mq/L© since bulges contribute typical 30% 
of the galaxy light in B (Schechter & Dressier 1987). Evidently, if e = 
0.1, we have already found about three times that value. This implies that 
either e is smaller than 0.1, or that quasars do not make up all of the AGN 
population. 

8. Are Supermassive Black Holes Ubiquitous? 

They certainly have not been found in every case that has been looked. Ko- 
rmendy has undertaken a systematic survey of a modest sample of galaxies 
(E-Sb), and his detection rate has been about 20% (KR). But, of course, 
many factors conspire against the detection of MDOs, and this estimate 
should be regarded as firm lower limit. If one takes seriously the M,-Mbui 
relation described above, it is possible that every bulge contains a SMBH 
with an appropriately scaled size. This view is supported by the statisti- 
cal analysis of Magorrian et al. (1998). In fact, the detection of an MDO 
in NGC 4945 (§ 4.2) and the presence of a bona fide AGN in NGC 4395 
indicate that perhaps even some galaxies without bulges may have SMBHs. 
Additional support for this picture comes from the growing evidence 
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that nonstellar nuclear activity is very common in galaxies, much more so 
than conventionally believed based on the statistics of bright AGNs and 
quasars. A recent spectroscopic survey of a large, statistically complete 
sample of nearby galaxies finds that over 40% of all bright {Bt < 12.5 mag) 
galaxies have nuclei that can be classified as "active," and the percentage 
is even higher among early- type systems (E-Sbc), approaching 50%-75% 
(Ho, Filippenko, & Sargent 1997). Most of the nearby AGNs have much 
lower luminosities than traditionally studied active galaxies, and a greater 
heterogeneity in spectral types is found (low-ionization nuclei, or LINERs, 
are common, for example), but the evidence is overwhelming that many of 
these nuclei are truly accretion-powered sources (see Filippenko 1996; Ho 
et al. 1997). Moreover, intrinsically weak, compact radio cores are known 
to be present in a significant fraction of elliptical and SO galaxies (Sadler, 
Jenkins, & Kotanyi 1989; Wrobel & Heeschen 1991), almost all of which 
spectroscopically qualify as AGNs (Ho 1998). 

Within the conventional AGN paradigm, the observed widespread nu- 
clear activity implies that SMBHs are a generic component of many, per- 
haps most, present-day bulge-dominated galaxies, consistent with the pic- 
ture emerging from the kinematic studies. This is a remarkable statement. 
It implies that SMBHs should not be regarded as "freaks of nature" that 
exist in only a handful of galaxies; rather, they must be accepted and under- 
stood as a normal component of galactic structure, one that arises naturally 
in the course of galaxy formation and evolution. 

9. Some Implications and Future Directions 

The SMBH hunting game is rapidly becoming a rather mature subject. 
I think we have progressed from the era of "the thrill of discovery" to a 
point where we are on the verge of using SMBHs as astrophysical tools. 
In this spirit, let me remark on a few of the ramifications of the existing 
observations and point out some of the more urgent directions that should 
be pursued. 

A. The M.-Mhiii relation. The apparent correlation between the mass of 
the central BH and the mass of the bulge, if borne out by future scrutiny, 
has significant implications (see below). Prom an observational point of 
view, the highest immediate priority is to populate the M.-M^ui diagram 
with objects spanning a wide range in luminosity, with the eventual aim of 
deriving a mass function for SMBHs. The samples should be chosen with 
the following questions in mind. (1) Is the apparent trend a true correlation 
or does it instead trace an upper envelope? (2) If the relation is real, is it 
linear? (3) What is the magnitude of the intrinsic scatter? And (4) is there 
a minimum bulge luminosity (mass) below which SMBHs do not exist? 
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In the near future, the most efficient way to obtain mass measure- 
ments for relatively large numbers of galaxies is to exploit the capabilities 
of STIS on HST. Several large programs are in progress. Although VLBI 
spectroscopy of H2O masers delivers much higher angular resolution, this 
technique is limited by the availability of suitably bright sources. Condi- 
tions which promote H2O megamaser emission evidently are realized in only 
a tiny fraction of galaxies (Braatz, Wilson, & Henkel 1996). 

B. The formation of SMBHs. The M,-Mbui relation offers some clues to 
the formation mechanism of SMBHs. How does a galaxy know how to ex- 
tract a constant, or at least a limiting, fraction of its bulge mass into a 
SMBH? An attractive possibility is by the normal dynamical evolution of 
the galaxy core itself. The spheroidal component of nearby galaxies can 
attain very high central stellar densities — up to 10^ Mq pc^'^ (Faber et 
al. 1997) — and some with distinct nuclei have even higher concentrations 
still (Lauer et al. 1995). Although most galaxy cores are unlikely to have 
experienced dynamical collapse (Kormendy 1988c), the innermost regions 
have much shorter relaxation times, especially when considering a realistic 
stellar mass spectrum because the segregation of the most massive stars 
toward the center greatly accelerates the dynamical evolution of the sys- 
tem. Lee (1995, and these proceedings) shows that, under conditions typical 
of galactic nuclei, core collapse and merging of stellar-size BHs can easily 
form a seed BH of moderate mass. Alternatively, the seed object may form 
via the catastrophic collapse of a relativistic cluster of compact remnants 
(Quinlan Sz Shapiro 1990). In either case, subsequent accretion of gas and 
stars will augment the central mass, and, over a Hubble time, may produce 
the distribution of masses observed. It is far too premature to tell whether 
SMBHs form through the secular evolution of galaxies, as suggested here, 
through processes associated with the initial formation of galaxies (e.g., 
Rees 1984, 1998; Silk & Rees 1998), or both, depending on the galaxy type 
(elliptical vs. spiral galaxies). But the stage is set for a serious discussion. 
More sophisticated modeling of the growth of SMBHs that take into ac- 
count a wider range of initial conditions in galactic nuclei (e.g., relaxing 
the adiabatic assumption or adopting more realistic density profiles for the 
stellar distribution) may eventually yield testable predictions (see, e.g., Sti- 
avelh 1998). 

C. Infi,uence of SMBHs on galactic structure. Norman, May, & Van Albada 
(1985) showed through A^-body simulations that a massive singularity in 
the center of a triaxial galaxy destroys the box-like stellar orbits and hence 
can erase the nonaxisymmetry, at least on small scales. This has several 
important consequences. First, it implies that the presence of a SMBH can 
influence the global structure and dynamics of galaxies. Second, the secular 
evolution of the axisymmetry of the central potential points to a natural 
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mechanism for galaxies to self-regulate the transfer of angular momentum 
of the gas from large to small scales. This negative-feedback process may 
limit the growth of the central BH and the accretion rate onto it, and 
hence may serve as a promising framework for understanding the physical 
evolution of AGNs. Merritt & Quinlan (1998) find that the timescale for 
effecting the transition from triaxiality to axisymmetry depends strongly on 
the fractional mass of the BH; the evolution occurs rapidly when M,/Mbui 
^ 2.5%, remarkably close to the observational upper limit (Fig. 8a). 

D. The oriQin of central cores. It is not understood why giant ellipticals 
have such shallow central light profiles. "Cores" do not develop naturally 
in popular scenarios of structure formation, and even if they form, they are 
difficult to maintain against the subsequent acquisition of the dense, central 
regions of satellite galaxies that get accreted (Faber et al. 1997). Moreover, 
the very presence of a SMBH, whether it grew adiabatically in a preexisting 
stellar system or the galaxy formed by violent relaxation around it (Stiavelli 
1998), ought to imprint a more sharply cusped light profile (see § 2) than 
is observed. An intriguing possibility is that cores were created as a result 
of mergers, where one or more of the galaxies contains a BH. As the single 
(Nakano & Makino 1998) or binary (Makino 1997; Quinlan &; Hernquist 
1997) BH sinks toward the center of the remnant due to dynamical friction, 
it heats the stars, thereby producing a "fluffy" core. If this interpretation is 
correct, it would provide a simple, powerful tool to diagnose the formation 
history of galaxies. 

E. Willi are Lhe black holes so black? It has been somewhat puzzling how 
the BHs can remain so dormant. No doubt the dwindled gas supply in the 
present epoch, especially in ellipticals, is largely responsible for the inac- 
tivity. Yet the accretion rate cannot be zero; even in the absence of inflow 
from the general interstellar medium, some gas is shed through normal 
mass loss from the innermost evolved stars, and occasionally such stars get 
tidally disrupted (see below). If SMBHs are indeed present, the radiative 
efficiency of the accretion flow must be orders of magnitude lower than that 
of "standard" optically thick, geometrically thin disks. Such a situation may 
be realized in accretion flows where advection becomes important when the 
accretion rate is highly sub-Eddington (Narayan &: Yi 1995; Abramowicz et 
al. 1995; Nakamura et al. 1996; Chakrabarti 1996). Sgr A* at the Galactic 
Center has a bolometric luminosity of only ~10^^ ergs s~^, or Lboi/^Edd 

3x10"^ (Narayan, Yi, & Mahadevan 1995); in the case of the LINER 
nucleus of M81 (Ho, Fihppenko, & Sargent 1996), Lboi ~ 10^^ ergs 
and i^boi/-^Edd ~ 10~^. The spectral energy distributions emitted by both 
objects differ dramatically from those of luminous AGNs and can be ap- 
proximately matched by advective-disk models. 

F. Tidal disruption of stars. The prevalence of SMBHs suggested by the 
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existing evidence predicts a relatively high incidence of tidal disruptions 
of stars as they scatter into nearly radial orbits whose pericenters pass 
within the tidal radius of the BH (Rccs 1998, and references therein). For a 
typical stellar density of 10^ stars pc'^, M, = 10^-10^ Mq, and a = 100- 
300 km s~^, a solar-type star will be disrupted once every 10^-10^ years. 
(BHs more massive than 10^ ]\I : will swallow the star whole.) Roughly 
half the debris becomes unbound and half gets captured into an accretion 
disk which undergoes a bright flare (~ lO^'^ Lq) lasting a few months to a 
year. The spectrum is expected to be mainly thermal and to peak in the 
extreme- UV and soft X-rays. The contribution to the near-UV and optical 
bands is uncertain; it depends on assumptions concerning the geometry of 
the accretion disk (thick or thin) and on whether an optically thick envelope 
can form. For plausible parameters, Ulmer (1998) estimates that a 10^ 
BH will produce a flare with an absolute magnitude of about -20 in U 
and -18.5 in V . The realization that SMBHs may be even more common 
than previously thought provides fresh motivation to search for such stellar 
flares; some observational strategics are mentioned by Rees (1998). Here, I 
wish to stress that quantifying the rate of stellar disruptions can be used as 
a tool to study the demography of SMBHs out to relatively large distances 
and hence should be regarded as complementary to the kinematic searches. 
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